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日 地 Halo 轨道 的 多 约束 转移 轨道 分 层 向 分 修正 设计 
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摘 要 : 针对 日 - 地 系统 LI 点 ( 简称 SEL1 点 ) Halo 轨道 转移 轨道 设计 中 存在 的 多 约束 与 初 值 敏感 性 问题 ， 
提出 一 种 基于 分 层 微分 修正 与 初 值 多 项 式 的 设计 方法 。 首 先 定义 平 动 点 转移 轨道 设计 过 程 中 存在 的 约束 条 件 , 然 
后 根据 不 同 的 终端 约束 条 件 ,重点 给 出 了 同时 考虑 轨道 高 度 、 轨 道 倾 角 、 升 交点 赤 经 与 航 迹 角 等 多 约束 条 件 下 的 分 
层 微分 修正 方法 。 通 过 分 析 约 束 变量 与 控制 变量 之 间 的 关系 ,得 到 能 够 解决 微分 修正 初 值 问题 的 初 值 表 达 式 。 最 
后 在 多 约束 条 件 下 设计 了 从 轨道 高 度 为 200 km 的 地 球 停泊 轨道 到 SEL1 点 Halo 轨道 的 转移 轨道 。 仿 真 结 果 表 
明 , 分 层 微 分 修正 方法 能 够 处 理 多 约束 问题 , 且 初 值 表 达 式 可 以 为 微分 修正 提供 良好 的 初始 条 件 , 从 而 保证 算法 收 
你 ,方法 具有 较 好 的 实用 性 。 
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Abstract: To cope with the orbit transfer problem with multiple constraints and sensitive initial values around the 
libration point L1 of the Sun 下 arth system (SEL] point), a transfer trajectory design method based on hierarchical 
differential corrections and initial value polynomial is proposed. At first, constraints are analyzed in the process of designing 
the transfer trajectory around the libration point. According to different terminal constraints, the hierarchical differential 
correction method is proposed by considering orbit altitude, orbit inclination, right ascension of ascending node and track 
angle simultaneously. Then, by analyzing relation between constraint variables and control variables, the initial value 
expression is got to solve the initial value problem of differential corrections. Finally, transfer trajectory from the earth 
parking orbit whose altitude is 200 km to Halo orbit around SELI1 point is designed under the condition of multi-eonstraints. 
The simulation results indicate the hierarchical differential correction method can deal with the problem with multiple 
constraints and initial value polynomial can provide appropriate initial conditions for differential corrections so as to 
guarantee convergence of the algorithm. In conclusion, this method has good practicality. 
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I 可 将 探测 器 定点 于 平 动 点 附近 ”。 例 如 , 当 航 
器 运行 在 Ll 点 附近 的 Halo 轨道 时 ,能 够 实现 对 日 
日 地 系统 Ll 点 ( 简称 SEL1 点 ) 是 位 于 日 地 之 间 0 并 可 为 探测 深 空 环境 提供 最 佳 的 落脚 点 。 因 
的 共 线 平 动 点 ,利用 SELl 点 相对 主 天 体 不 变 的 几何 。 此 ,在 考虑 时 间 与 燃 耗 等 条 件 下 ,探测 器 如 何 顺利 地 
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到 达 平 动 点 附近 轨道 成 为 了 现今 深 空 领域 研究 的 热 
点 问题 。 目 前 ,已 有 多 个 探测 器 利用 SEL1 点 Halo 轨 
道 执 行 了 科学 探测 任务 。 具 体 地 ,分 别 于 1978 年 和 
1994 年 发 射 的 ISEE3 及 Wind 探测 器 ,通过 SELL 点 
附近 的 Halo 轨道 ,观测 了 太阳 活动 变化 对 近 地 空 间 
环境 的 影响 ”。 于 1995 年 发 射 的 SOHO 观测 台 , 运 
行 在 Halo 轨道 上 观测 分 析 了 太阳 内 部 结构 变化 , 进 
一 步 研 究 了 太阳 风 的 形成 与 加 速 过 程 ”。 而 NASA 
研制 的 Genesis 探测 器 同样 利用 SEL1 点 Halo 轨道 完 
成 了 采集 返回 太阳 风 样 本 任务 ” 。 

针对 平 动 点 转移 轨道 设计 ,文献 5] 详细 地 分 
析 了 平 动 点 附近 同 、 异 宿 连接 的 力学 机 制 并 论述 了 
平 动 点 理论 在 Halo 轨道 的 转移 轨道 设计 中 的 应 用 。 
Zazzera 等 结合 Poincare 截面 与 不 变 流 形 设 计 了 无 
约束 条 件 下 的 平 动 点 转移 轨道 8 。 李 明 涛 等 以 小 
振幅 Halo 轨道 为 目标 轨道 并 利用 不 变 流 形 分 别 对 
不 合约 束 条 件 与 考虑 轨道 高 度 、 轨 道 倾角 或 升 交点 
赤 经 时 的 L1 点 转移 轨道 进行 了 分 析 ”“” 。 文 献 

日 将 序 优化 理论 与 微分 修正 算法 结合 ,在 仅 考虑 

地 球 停泊 轨道 高 度 约束 条 件 下 , 设计 了 SELL 点 
Halo 轨道 的 两 脉冲 转移 方案 。Qi " 与 雷 汉 伦 等 
等 结合 不 变 流 形 与 小 推力 技术 ,利用 混合 优化 算法 ， 
分 别 设计 了 满足 状态 矢量 约束 的 地 -月 与 日 - 火 系 
平 动 点 转移 轨道 。Howell 等 &2 利用 两 级 微分 修正 
法 设计 了 特定 轨道 倾角 及 升 交 点 赤 经 约束 下 的 
Genesis 返回 任务 轨道 。 但 两 级 微分 修正 存在 分 段 
优化 计算 时 间 长 与 初 值 精度 要 求 高 的 问题 。 

综 上 所 述 ,以 往 研 究 常常 局 限于 无 约束 或 两 种 
以 下 约束 ( 简称 少 约 束 ) 的 简化 模型 , 基于 简化 模型 
必 闻 增 大 轨道 设计 误差 ,难以 满足 越 来 越 复 杂 的 任 
务 要 求 。 并 且 转 移 轨 道 末 端 少 约束 条 件 无 法 完全 确 
定 地 球 停泊 轨道 根 数 ,在 工程 应 用 中 存在 着 局 限 性 。 
与 此 同时 ,在 轨道 设计 过 程 中 ,以 往 研 究 并 未 给 出 详 
细 的 初 值 选取 依据 。 

本 文 针对 以 上 问题 ,提出 了 相应 解决 方法 。 首 
先 , 在 轨道 设计 过 程 中 同时 考虑 轨道 高 度 、 轨 道 倾 
角 、 升 交点 赤 经 及 航 迹 角 等 多 种 约束 条 件 ( 简称 多 
约束 ) , 从 而 能 够 准确 地 确定 圆 形 地 球 停泊 轨道 。 
在 求解 多 约束 情况 下 的 转移 轨道 设计 问题 时 , 提出 
了 分 层 微分 修正 搜索 算法 ,解决 了 一 次 修正 算法 容 
易 出 现 不 收敛 的 问题 , 并且 通 过 分 析 状 态 转 移 矩 阵 


在 微分 修正 中 的 作用 ,在 算法 中 添加 上 比例 因子 加 快 
收敛 ,提高 了 计算 效率 。 其 次 ,为 了 解决 微分 修正 算 
法 的 初 值 敏感 性 问题 ,研究 了 目标 轨道 入 轨 点 不 同 
机 动 速度 增 量 大 小 及 方向 对 转移 轨道 设计 的 影响 ， 
给 出 了 初 值 多 项 式 , 保 证 修正 算法 收敛 。 通 过 有 具体 
算 例 分 析 , 有 效 地 设计 了 从 轨道 高 度 为 200 km 的 地 
球 停泊 轨道 出 发 到 SELI 点 Halo 轨道 不 同 入 轨 点 的 
转移 轨道 ,对 深 空 探测 器 的 轨道 设计 具有 借鉴 意义 。 


1 圆 型 限制 性 三 体 模型 及 设计 变量 


本 文 以 日 `. 地 及 航天 器 构成 的 圆 型 限制 性 三 体 
模型 进行 分 析 , 研究 航天 器 在 两 个 主 天 体 引 力作 用 


下 的 运动 状态 。 根 据 能 量 法 并 对 参数 进行 无 量 纲 化 
后 ,航天 器 在 会 合 坐标 系 下 的 动力 学 方程 为 ": 
xX-2y=-U, 
Yt (1) 
z=-U, 
式 中 
机 业 > 1] - 
OS 
ni (x+ 内 + +2 
Fyn (x -1+H ty + 


(2) 
U,，U,,U, 分 别 表示 有 效 势能 U 关于 *,y;z 的 偏 导 
数 。 
为 了 清晰 地 描述 转移 轨道 与 初 值 表达 式 的 设计 
过 程 ,首先 定义 转移 轨道 设计 时 使 用 的 控制 变量 与 
约束 变量 ,如 图 1 所 示 。 
一 一 ~ 转移 轨道 


i A 


1 转移 轨道 设计 变量 
Fig.1 Definition of transfer trajectory design variables 
为 了 充分 利用 目标 Halo 轨道 的 信息 并 保证 转 
移 轨道 末端 满足 约束 条 件 , 本 文采 用 反 向 积分 方法 
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求解 转移 轨道 ,即将 Halo 轨道 作为 转移 轨道 的 初始 
点 ,而 终点 位 于 地 球 停泊 轨道 上 。 因 此 ,在 微分 修正 
过 程 中 ,控制 变量 选取 Halo 轨道 入 轨 点 ( 简称 HOI 
点 ) 处 三 轴 速 度 增 量 刀 = 从 V,,AV,,AV」“" 及 飞行 
时 间 7 。 在 停泊 轨道 逃逸 点 ( 简称 TTI 点 ) 处 ,为 了 
直观 地 描述 地 球 停泊 轨道 参数 ,选取 约束 变量 处 ,= 

有 ,yi, 人 ,分 别 为 轨道 高 度 、 航 迹 角 、 轨 道 倾角 
与 升 交点 赤 经 。 

同时 ,为 了 保证 控制 变量 具有 良好 的 初 值 ,从 而 
使 得 微分 修正 算法 收敛 ,对 初 值 多 项 式 设计 过 程 中 
使 用 到 的 变量 进行 定义 ( 见 图 1) 。 在 HOI 点 处 ,将 
速度 增 量 进行 投影 ,其 中 AV, ,AV., 分 别 为 xy 平面 
与 z 轴 上 的 速度 增 量 , a,6 为 各 坐标 平面 内 由 速度 
增 量 构 成 的 角度 。 在 TTI 点 处 , y 为 相对 地 球 的 飞 
行 航 迹 角 , 6 为 位 置 矢量 在 xy 平 面 内 投影 ,与 * 轴 
的 夹 角 。 始 末端 点 变量 的 数值 可 分 别 通过 初 值 点 与 
终点 的 状态 矢量 进行 求解 。 

由 于 地 球 停泊 轨道 升 交 点 赤 经 与 轨道 倾角 
定义 在 地 心 赤道 惯性 系 下 ,而 求解 航天 器 的 动力 
学 方程 位 于 会 合 坐标 系 下 ,因此 需要 将 不 同 坐标 
系 下 的 状态 量 转换 到 同一 坐标 系 中 。 本 文 将 会 
合 系 下 的 状态 矢量 转换 到 地 心 赤道 惯性 系 下 , 对 
轨道 倾角 与 升 交点 赤 经 进行 求解 ,首先 定义 图 2 
所 示 坐 标 系 。 


2 坐标 系 定义 


Fig.2 The definition of coordinate system 


SX,Y.Z, 为 日 心 黄道 坐标 系 ; Bxyz 为 质心 会 合 
坐标 系 ; Ex.y,z. 为 地 心 会 合 坐标 系 。 同 时 , 假设 地 
球 相 对 于 太阳 的 升 交 点 赤 经 为 2, 、 轨 道 倾角 为 i, 及 
纬度 幅 角 为 9 , 地 球 轨道 根 数 取 值 通过 星 历 获 得 。 
则 地 心 黄道 坐标 系 与 地 心 会 合 系 之 间 的 坐标 转换 矩 
阵 为 


coco — SpCisg, — CoSg, 一 SCiCco 5330 
C= Sp.Cg, + CoCisg, — sosg, tcocice, -Sico, 
S750, Si Cg Ci, 


(3) 
式 中 : cs 分 别 表示 余弦 值 与 正弦 值 。 对 式 (2) 关 
于 时 间 t 求 导 可 得 速度 矢量 表达 式 。 其 次 ,假设 黄 
道 面 与 地 球 赤道 面 之 间 的 夹 角 为 。 。 最 终 , 状态 量 
由 质心 会 合 坐 标 系 转换 到 地 心 赤 道 惯性 系 为 


x a 
多 y 
入 四 R pa 
x 0,s Rix 
y. 7 
i . 
(4 -DC ] 
cp -1) Cr +s(l -NA)C 
se( -1) Cr +e(p -1)C 
(1 -1) Co 
co( —1)C» +s(l -nh) Cy 
LHC +ces(p -1) C1 
(4) 
式 中 : 
Ci Gi Cis 
R=|cC -SC CaC»y -soC3 CC — ssC3 
Ss Ca +L SC +e Gs 8 C0 +e Ca 


CG; 为 坐标 转换 矩阵 中 第 i 行 、 第 j 列 元 素 , 则 轨道 倾 
角 与 升 交点 赤 经 的 表达 式 满 足 


Ney)e YeXe 


i = arceos| 


(2s a) + (20, = Weg 


2 


+ (xy。 — ye) (5) 


NZe ZeNe 


(2 = arceos| 


(re ma 
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应 注意 ,在 求解 升 交 点 赤 经 时 , 取 值 存在 象限 问题 ， 


需要 分 类 讨论 。 
2 约束 条 件 下 的 转移 轨道 设计 方法 


2.1 末端 约束 条 件 分 析 

转移 轨道 未 端 常常 施加 约束 条 件 以 满足 实际 任 
务 的 要 求 。 同 时 ,针对 圆 形 停泊 轨道 ,轨道 根 数 包 含 
轨道 高 度 、 轨 道 倾 角 及 升 交点 赤 经 。 因 此 ,为 了 准确 
地 确定 地 球 停泊 轨道 ,本 文 同时 选取 轨道 高 度 、 轨 道 
倾角 、 升 交点 赤 经 及 航 迹 角 作 为 约束 条 件 。 结 合 
述 四 种 约束 条 件 ,求解 控制 变量 与 约束 变量 关系 式 
的 步骤 如 下 : 

1) 根据 转移 轨道 设计 要 求 , 选 定 控制 变量 。 通 
过 分 析 , 本文 选取 控制 变量 为 

X, = WAV,,AV,,AV.,N" (6) 

2) 针对 轨道 未 端 施加 的 约束 条 件 , 利 用 末端 状 
态 量 计算 约束 条 件 的 实际 值 ,并 与 期 望 值 X, 进行 比 
较 得 到 约束 矢量 , 即 
h, 加 h 
Yr Ya (7) 
Ua 
-0 

3) 通过 分 析 控 制 变量 与 约束 变量 的 关系 ,计算 
约束 条 件 相 对 于 控制 变量 的 偏 导 数 , 即 雅克 比 和 矩阵 ， 
可 表示 为 


-ogAV, 9AV, 9AV. 97] 

4) 在 实际 求解 时 ,常常 以 航 迹 角 为 雳 作为 积 
分 终止 条 件 ,能 够 简化 计算 并 且 得 到 合适 的 转移 
轨道 , 即 地 球 停泊 轨道 半径 不 小 于 地 球 半径 。 
此 , 设 定 6y = yj -yv =0。 同 时 ,结合 雅克 比 矩 
阵 可 推导 出 函数 67 = 67( 6AV,,6AV,,6AV.) , 消 
去 时 间 变 量 后 得 到 控制 变量 与 约束 变量 关系 式 
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为 
6h, An 42 4A]1T6AV, 
6i |= |4, 42 2 | 6AV (9) 
680 LA, A A;JL8AV 
以 轨道 高 度 为 例 
oh 
1 oAV, 6y oAV, 
oT 
ah, 
oh, aT 
hr sy (10) 
aAV, ay 9AV, 
oT 
ohy 
沂 E oh _o7 97 
8 oAV. 9y 9AF 
oT 


轨道 倾角 与 升 交点 赤 经 约束 同样 具有 类 似 表达 
式 。 需 要 注意 ,在 具体 地 求解 过 程 中 , 应 考虑 会 合 4 
标 系 与 地 心 赤 道 惯性 系 之 间 的 坐标 转换 矩阵 。 
2.2 微分 修正 初 值 多 项 式 

微分 修正 作为 一 种 局 部 迭代 方法 ,算法 能 否 收 
敛 到 期 望 值 ,很 大 程度 上 取决 于 初 值 是 否 准确 ,而 在 
以 往 的 研究 中 并 没有 给 出 确切 的 初 值 及 选取 依据 。 
本 文通 过 分 析 初 始 控制 变量 与 未 端 约束 变量 之 间 的 
关系 , 对 数据 进行 多 项 式 拟 合 ,最 终 得 到 包含 轨道 倾 
角 与 飞行 时 间 等 参量 的 初 值 表达 式 , 为 轨道 设计 时 
确定 初 值 提供 方便 。 

具体 地 ,设计 过 程 中 限定 HOI 处 的 法 向 速度 增 
量 大 小 , 即 A『 € 广 30 m/s,30 m/s| 并 以 2.5 m/s 
的 增 量 增加 ,针对 特定 法 向 速度 增 量 大 小 , 通过 改变 
xy 平面 的 速度 增 量 与 6 角 取 值 , 从 而 得 到 满足 转移 
轨道 末端 位 于 地 球 附 近 等 要 求 的 大 量 数据 。 基 于 数 
据 并 分 析 HOI 点 处 不 同 速度 增 量 大 小 与 飞行 时 间 、 
轨道 倾角 等 参数 之 间 的 关系 ,并 结合 图 1 中 定义 的 
角度 ,对 数据 进行 拟 合 后 ,最 终 得 到 包含 轨道 根 数 与 
飞行 时 间 等 参量 的 初 值 表达 式 。 

3 给 出 了 B 角 与 飞行 时 间 之 间 的 关系 。 由 图 
可 知 , 在 不 同 的 法 向 速度 增 量 大 小 条 件 下 ,曲线 变化 
趋势 一 致 并 近似 符合 三 次 函数 的 特性 。 因 此 , 基于 
三 次 函数 进行 构造 , 6 角 与 飞行 时 间 之 间 满 足 

了 = 99. 486 + 0.5198 + 0.01498” - 0. 00028° 
(11) 
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飞行 时 间 / 大 


a0 ”和 人 
BA) 


3 B 角 与 飞行 时 间 变 化 曲线 


Fig.3 The relation curve between B angle and time of flight 


10 0 30 
P/O) 


4 ”速度 增 量 与 6 角 之 间 的 关系 


Fig.4 The relation between B angle and the velocity increment 


4 描述 了 速度 增 量 治平 面 及 轴线 分 量 大 小 与 B 
角 之 间 的 关系 曲面 ,其 中 的 点 表示 仿真 数据 ,具有 二 次 
曲面 特性 ,针对 此 特性 可 得 到 式 ( 12) 所 示 的 表达 式 
AV,, = -3.958B + 0. 1348” + 0.0178BAV. - 
0. 190AV. + 0.00604AV + 61.553 (12) 
由 式 ( 11) 与 式 ( 12) 可 知 , 在 求解 过 程 中 需要 确 
定 法 向 速度 增 量 大 小 AF 与 a 角 的 取 值 ,而 不 同 的 
飞行 时 间 及 转移 轨道 末端 约束 条 件 对 取 值 存在 着 一 
定 的 影响 。 因 此 , 基于 多 变量 对 速度 增 量 大 小 AV. 
与 a 角 的 表达 式 进 行 拟 合 求解 ,满足 
i = -0.1758 + 0.002668”- 0. 55242 + 0. 00988(2 - 
0. 0218B80 + 0. 000193B(2 - 5.721 x 
10 “(7 + 33. 02 
AT = (4.6537 - 0.01537 + 0.0183. - 
0. 000188(2 - 311. 646 - B) /0. 0122 
a = 23.007 -0.5847 - 0.00164i - 0.00362AV. - 
0.08158 + 0.003627T - 7.80 x 1072 ~ 0.003058” 
( 13) 


式 中 :角度 取 值 范围 为 8 s 0D°,30°] ,2 < [-28°， 
60°] ,此 处 升 交点 赤 经 选取 负 值 是 为 了 保证 得 到 良 
好 的 初 值 表达 式 。 

表 1 非 线 性 拟 合 相 关 参 数 


Table 1 Parameters related to nonlinear fitting 

拟 合 画 数 相关 系数 平方 RMSE 

i( Q,B) 0.9956 0.6394 

AV,,(B,AV.) 0.9982 0.4572 

AV.( 7,0,8) 0. 9955 0.3418 
a( TOF, AV. ,i,B) 0.9917 2.09 x10-8 


针对 微分 修正 初 值 多 项 式 的 设计 结果 ,本 文 以 
最 大 残 差 模 、 相 关系 数 平方 值 ” 与 均 方 根 误 差 
RMSE 为 依据 判断 拟 合 程度 的 好 坏 。 其 中 , 图 3 拟 
合 曲线 对 应 的 最 大 残 差 模 为 2. 9741 , 而 其 他 多 变量 
拟 合 画 数 的 相关 参数 如 表 1 所 示 。 通 过 非 线性 拟 合 
参数 与 仿真 校 验 可 知 , 微 分 修正 初 值 多 项 式 拟 合 效 
果 较 好 ,能 够 为 控制 变量 提供 良好 初始 条 件 并 保证 
算法 收敛 。 

2.3 ”分 层 微 分 修正 搜索 算法 

微分 修正 法 为 一 种 牛顿 迭代 逼近 期 望 值 算法 ， 
通过 不 断 修正 控制 变量 使 得 约束 变量 最 终 收 敛 于 期 
望 值 ” 。 当 同时 考虑 轨道 高 度 .轨道 倾角 、 升 交点 
赤 经 及 航 迹 角 多 个 约束 时 ,由 于 转移 轨道 末端 满足 
约束 条 件 的 程度 各 不 相同 ,相互 之 间 存 在 影响 ,最 终 
使 得 一 次 修正 算法 将 出 现 不 易 收 敛 的 问题 。 再 者 ， 
利用 遗传 算法 或 序列 二 次 规划 等 优化 算法 对 多 约束 
问题 进行 求解 时 ,存在 计算 时 间 长 等 问题 。 因 此 , 针 
对 以 上 间 题 ,结合 转移 轨道 末端 不 同 的 约束 条 件 , 利 
用 分 层 修正 算法 搜索 精确 轨道 ,使 得 转移 轨道 末端 
同时 满足 轨道 高 度 、 轨 道 倾角 、 升 交点 赤 经 及 航 迹 角 
的 多 种 约束 条 件 。 

第 一 层 修正 : 使 转移 轨道 未 端 靠 近 地 球 。 

在 初步 搜索 过 程 中 ,只 考虑 轨道 高 度 约束 。 由 
于 选取 的 停泊 轨道 为 近 地 轨 道 且 初次 积分 得 到 的 转 
移 轨道 末端 距离 地 球 较 远 ,通过 第 一 层 修正 能 够 保 
证 转移 轨道 末端 位 于 停泊 轨道 附近 ,为 施加 倾角 和 与 
升 交 点 赤 经 约束 委 定 基础 。 具 体 地 , 首先 对 初始 状 
态 量 进行 积分 ,判断 约束 条 件 Ah, 是 否 小 于 给 定 误 
差 , 若 满足 则 进行 下 一 层 微分 修正 。 否 则 ,微分 修正 
控制 变量 X = ”AV, ,AV,,AV]" ,直至 满足 轨道 末 
端 位 置 约束 条 件 。 

第 二 层 修正 : 考虑 轨道 高 度 与 倾角 约束 。 
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经 过 第 一 层 修正 后 ,保证 了 转移 轨道 末端 位 于 
地 球 附 近 , 此 阶段 进一步 细 化 修正 模型 ,同时 考虑 轨 
道 高 度 与 轨道 倾角 约束 。 以 第 一 层 满足 高 度 约束 的 
状态 量 为 初 值 进 行 积分 ,判断 轨道 未 端 高 度 与 倾角 
约束 Ah, ,Aii 是 否 小 于 误差 值 , 若 满足 则 进行 第 三 
层 修正 ,否则 进一步 对 控制 变量 循环 迭代 求解 ,使 得 
最 终 约束 收敛 于 期 望 值 。 

第 三 层 修 正 : 求解 精确 SEL1 点 Halo 轨道 转移 
轨道 。 

通过 前 两 层 微分 修正 ,转移 轨道 末端 高 度 与 倾 
角 已 经 接近 期 望 值 , 但 仍 不 满足 升 交点 赤 经 约束 。 
因此 ,第 三 层 修正 同时 考虑 轨道 高 度 、 轨 道 倾 角 及 升 
交点 赤 经 约束 , 最 终 获得 满足 多 约束 条 件 的 精确 
SEL1 点 Halo 轨道 转移 轨道 。 同 样 地 , 对 上 一 层 状 
态 量 进行 积分 ,迭代 修正 控制 变量 使 得 轨道 末端 同 
时 满足 轨道 高 度 、 倾 角 及 升 交 点 赤 经 多 约束 条 件 。 
综 上 所 述 , 分 层 微分 修正 流程 如 图 5 所 示 。 

另 一 方面 ,为 了 保证 算法 具有 较 好 的 收敛 效果 
且 缩 短 计算 时 间 ,本文 对 传统 微分 修正 算法 中 迭代 
修正 量 与 状态 转移 矩阵 分 别 进行 改进 。 

首先 针对 迭代 修正 量 ,为 了 保证 算法 收敛 ,对 每 
一 步 的 迭代 修正 量 乘 以 式 (14) 所 示 的 步 长 因子 。 
主要 功能 为 : 由 于 微分 修正 初始 阶段 ,实际 值 与 期 望 
值 差距 较 大 ,通过 步 长 因子 能 够 保证 较 小 的 修正 量 ， 
防止 一 步 “器 过 ”期 望 值 而 导致 算法 发 散 。 同 时 , 随 
着 迭代 修正 次 数 不 断 增加 , 步 长 因子 逐步 增 大 并 最 
终 趋 近 于 1 ,不 会 影响 收敛 速度 。 即 

fa (14) 

式 中 : a,5 为 常数 且 取 值 a = 1,5 = 0. 85 ,x 在 每 步 
迭代 修正 中 以 w = 0. 1 增 量 增加 。 

而 针对 状态 转移 矩阵 ,通过 仿真 校 验 , 在 微分 修 
正 过 程 中 , 对 状态 转移 矩阵 进行 比例 变换 后 ,微分 修 
正 迭 代 次 数 将 会 减少 。 有 具 体 地 ,上 比例 系数 可 根据 当 
前 修正 量 与 前 一 时 刻 修正 量 进行 比较 后 增加 或 减 
小 。 以 只 考虑 轨道 高 度 约束 的 第 一 层 修正 为 例 , 定 
义 判 别 依据 5AA = Ah, - Ah,_， ; 则 比例 系数 及 可 表 


B= -gqg, 6Ah<s0,i>1 (15) 
B= rg 6Ah >0,:>1 


式 中 : 9 为 取 值 小 于 0.1 的 实数 ,本 文中 选取 gq，= 
0. 025 ,而 ; 为 微分 修正 迭代 次 数 。 第 二 层 与 第 三 层 
修正 分 别 选取 轨道 倾角 与 升 交 点 赤 经 为 判别 依据 。 


根据 初 们 表达 式 确定 HOI 
点 满足 任务 要 求 的 速度 


增 量 AV,,AV,a 


AV', AV ,AV. 


团 型 限制 性 王 体 模型 下 
轨道 积分 


连 代 修正 控制 变 世 


于 型 限制 性 下 体 模 型 下 
轨道 积分 


迁 代 修 由 控制 变量 
AV,AV,,AV. 


Ah, < tol, 
Ai < tol; 


司 型 限制 性 一 体 模 型 下 


机 道 积分 


迭代 修 下 控制 变 昌 
AV,AV,,AV, 


图 5 转移 轨道 分 层 修正 流程 图 
Fig.5 Flow chart of hierarchical differential corrections 


of transfer trajectory 


3 仿真 结果 分 析 


3.1 多 约束 条 件 的 转移 轨道 与 算法 比较 

本 文 以 ISEE -3,SOHO 等 太阳 观测 器 的 目标 轨 
道 为 研究 对 象 ,即日 地 系统 Ll 点 附近 法 向 幅 值 为 
120000 km 的 Halo 轨道 ,利用 微分 修正 初 值 多 项 式 
为 控制 变量 提供 初 值 ,设计 了 多 约束 条 件 下 的 平 动 
点 转移 轨道 ,并 针对 转移 轨道 末端 不 同 约束 条 件 下 
的 飞行 时 间 与 速度 增 量 进行 讨论 。 

在 求解 会 合 坐标 系 与 地 心 赤道 惯性 坐标 系 间 的 
坐标 转换 矩阵 时 , 选取 发 射 时 刻 为 2015 年 8 月 10 
日 ,对 应 的 地 球 轨道 根 数 满足 : 升 交点 赤 经 2，= 
347.95° ,轨道 倾角 i。 = - 0.00198° ,纬度 幅 角 
0. = - 31.08°。 地 球 停泊 轨道 要 素 及 任务 约束 要 
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求 如 表 2 所 示 。 其 中 , 升 交点 赤 经 在 式 ( 13) 给 定 角 
度 范围 内 取 值 ,而 轨道 倾角 的 选取 符合 我 国 发 射 场 
对 倾角 的 要 求 。 

表 2 任务 约束 条 件 


Table 2 The mission constraints 


参数 约束 条 件 
轨道 高 度 h /km 200 
偏心 率 e 0 
轨道 倾角 i /(“) 50 
升 交 点 赤 经 2 /(°) -10 
航 迹 角 y /(°) 0 
角度 B /(。) 20 
位 置 误差 110°7 
角度 误差 1x1074 


针对 以 上 约束 条 件 ,采用 分 层 微分 修正 方法 
进行 转移 轨道 设计 。 首 先 利用 第 一 层 微分 修正 
将 转移 轨道 未 端 轨道 高 度 限定 在 400 km。 然 后 
在 第 二 层 微分 修正 中 考虑 轨道 倾角 约束 ,使 得 未 
端 满足 轨道 高 度 300 km, 轨道 倾角 i = 40。 约 束 
条 件 。 最 后 通过 第 三 层 修正 得 到 同时 满足 轨道 
高 度 、 轨 道 倾 角 、 升 交点 赤 经 及 航 迹 角 多 个 约束 
条 件 的 精确 转移 轨道 。 
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6 考虑 轨道 高 度 约束 的 转移 轨道 


Fig.6 Transfer trajectory design with altitude constraints 


6 为 第 一 层 微分 修正 后 转移 轨道 在 会 合 坐标 
系 下 的 三 维 视 图 及 其 投影 图 。 在 此 层 修 正中 , 转移 
轨道 飞行 时 间 为 105. 5 天 , 地 球 处 机 动 速 度 增 量 大 
小 为 3208. 89 m/s, Halo 轨道 入 轨 机 动 速 度 增 量 大 
小 为 35.45 m/s。 

7 为 考虑 轨道 高 度 与 轨道 倾角 约束 后 的 转移 
轨道 。 施 加 轨道 倾角 约束 后 ,需要 消耗 一 定 的 速度 


增 量 值 对 第 一 层 转 移 轨道 进行 调整 。 此 时 , 转移 轨 
道 飞 行 时 间 为 103. 2 天 ,两 次 机 动 速度 增 量 大 小 分 
别 为 3233. 61 m/s 与 51.39 m/s。 


1.0 
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7 ”施加 轨道 倾角 后 的 转移 轨道 投影 


Fig.7 The projection of transfer trajectory with 
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8 ”精确 转移 轨道 


Fig.8 The precise transfer trajectory 


8 描述 了 同时 满足 轨道 高 度 、 轨 道 倾 角 、 升 交 
点 赤 经 及 航 迹 角 多 个 约束 的 精确 转移 轨道 。 首 先 ， 
探测 器 施加 3258. 23 m/s 大 小 的 速度 增 量 从 
200 km 高 度 地 球 停泊 轨道 逃逸 后 ,经 过 109.9 天 飞 
行 到 达 SEL1 点 附近 Halo 轨道 入 轨 点 ,再 次 施加 速 
度 垃 量 56. 09 m/s 进入 目标 Halo 轨道 。 

同时 ,为 了 比较 单 层 微分 修正 算法 与 分 层 微 
分 修正 算法 在 处 理 多 约束 条 件 时 存在 的 差异 。 
本 文 在 不 同 的 初始 条 件 与 约束 条 件 下 ,基于 初 值 
多 项 式 , 分 别 利 用 两 种 修正 算法 进行 仿真 ,并 分 
析 两 种 算法 在 多 约束 条 件 时 的 收敛 情况 与 迭代 
次 数 。 

本 文 针对 HOI 处 不 同 法 向 速度 增 量 大 小 , 即 
AV. =-5 m/s,0 m/s,5 m/s 三 种 情况 进行 分 析 , 约 
束 条 件 满足 : 轨道 高 度 h = 200 km ; 轨道 倾角 i = 
35° 并 以 3° 增 量 增加 直至 i = 50°; 升 交点 赤 经 
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人 2 = -10° 并 以 4° 垃 量 增加 直至 (2 = 10。。 因 此 ， 
每 个 法 向 速度 增 量 大 小 下 需要 分 析 36 种 不 同 的 约 
束 条 件 。 
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9 不同 约 束 条 件 下 单 层 ( 左 ) 与 分 层 ( 右 ) 
微分 修正 的 收敛 情况 
Fig.9 The astringency of monolayer ( left) and hierarchical 


(right) differential corrections with different constraints 


目 单 层 微分 修正 。 上 分 层 微分 修正 


收敛 算 例 数目 总 和 


不 同 取 值 A 


10 不 同 法 向 速度 条 件 下 ,收敛 算 例 数目 比较 图 
Fig. 10 The comparison of the number of convergent examples 


with different normal velocity 


9 与 图 10 分 别 给 出 了 三 种 法 向 速度 取 值 条 
件 下 ,利用 单 层 与 分 层 微 分 修正 算法 处 理 多 约束 问 
题 时 的 收敛 情况 与 对 应 的 收敛 算 例 数 目 ,图 9 中 点 
代表 满足 约束 时 对 应 的 轨道 倾角 与 升 交点 赤 经 。 同 
时 ,为 了 比较 两 种 修正 算法 的 计算 时 间 ,本文 以 友 代 
次 数 为 判断 依据 进行 分 析 。 如 图 11 所 示 , 描 述 了 法 
向 幅 值 A = - 5 m/s 时 ,两 种 算法 在 处 理 不 同 约束 
条 件 的 迭代 次 数 分 布 图 ,不 同 约束 对 应 的 迭代 次 数 
以 点 表示 ,其 中 编号 1 对 应 于 i = 35°,(2 = - 109) 
号 7 则 对 应 于 约束 i = 38",O = - 10° ,其 他 编号 以 
此 类 推 。 

由 图 9 与 图 10 可 知 , 在 某 些 约束 条 件 下 , 单 
层 修 正 算法 无 法 使 得 转移 轨道 末端 满足 给 定 的 


约束 条 件 , 而 分 层 微 分 修正 却 能 够 很 好 地 保证 算 
法 收敛 于 期 望 约束 值 。 同 时 ,从 图 11 可 见 ,在 满 
足 约束 的 条 件 下 , 单 层 与 分 层 微分 修正 算法 的 达 
代 次 数 相 近 , 但 在 部 分 情况 下 ,分 层 微分 修正 算 
法 的 返 代 次 数 优 于 单 层 微分 修正 迭代 次 数 ,其 主 
要 原因 是 ,在 某 些 条 件 下 ,为 了 保证 单 层 修正 算 
法 的 收敛 性 ,需要 式 ( 14) 中 取 值 较 小 ,进而 堆 
加 了 迭代 次 数 。 因 此 ,与 单 层 修正 算法 相 比 ,分 
层 修 正 算法 能 够 更 好 地 处 理 转移 轨道 设计 中 的 
多 约束 问题 。 


一 9 一 分 层 微分 修正 ~ 全 ~ 单 屋 微 分 修正 


33537 WHISNSTNI921232352729.3130.35 
算 例 编号 
11 法 向 幅 值 AV. = - 5 m/s ,迭代 次 数 比较 图 
Fig.11 The comparison of iterations with the normal 


amplitude AV, = -SS m/s 


3.2 不 同 约束 条 件 对 转移 轨道 的 影响 

为 了 验证 初 值 表达 式 在 微分 修正 过 程 中 的 有 效 
性 ,并 分 析 不 同 约束 条 件 取 值 对 转移 轨道 设计 的 影 
响 , 分 别 对 少 约束 ` 固 定 升 交 点 赤 经 与 固定 轨道 倾角 
的 多 约束 条 件 下 转移 轨道 进行 设计 分 析 。 同 时 , 结 
合 我 国 发 射 场 因 素 选 定 地 球 停泊 轨道 参数 , 即 轨道 
倾角 大 于 发 射 场 纬 度 、 东 南方 向 发 射 。 各 项 设计 参 
数 如 表 3、 表 4 及 表 5 所 示 。 

表 3 少 约束 条 件 下 的 转移 轨道 
Table 3 Transfer trajectory design with few constraints 


TTI 点 机 动 速度 HOI 点 机 动 速度 


增 量 /(m。s-!) 增 量 /(m。s-0) 


轨道 倾角 /(°) ”飞行 时 间 / 天 


30 106.8 3258.25 34.61 
40 108.9 3258.37 43.05 
50 107.3 3258.53 54.28 
60 107.4 3258.71 67.59 
70 108.3 3258.96 81.35 
80 109.7 3259.26 96.71 
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表 4 ”固定 升 交 点 赤 经 的 多 约束 转移 轨道 
Table 4 Transfer trajectory with multiple constraints 


for fixed RAAN 


升 交 点 ”轨道 倾 飞行 。 ”TTL 点 机 动 速度 HOI 点 机 动 速度 
赤 经 /(°) 角 /(*) 时间 /天 ” 增 量 /(m，s-!) 雹 量 /(m vs -0) 


30 110.2 3258.20 35.77 
40 100.9 3258.25 54.34 
50 109.9 3258.23 56.09 
让 60 114.8 3258.25 77.47 
70 118.1 3258.38 98.69 
80 120.7 3258.31 119.01 


表 5 固定 轨道 倾角 的 多 约束 转移 轨道 
Table 5 Transfer trajectory with multiple constraints for 


fixed inclination constraints 


轨道 倾 升 交 点 。 ”飞行 时 TTI 点 机 动 速度 HOI 点 机 动 速度 


角 /(°) 赤 经 /(*) ” 间 / 天 ” 增 量 /(m*s-!') 增 量 /(m vs- ) 
-10 110.2 3258.20 35.77 
0 104.7 3258.27 38.10 
10 114.8 3258.37 43.78 
20 110.1 3258.24 38.55 
30 116.1 3258. 84 73.09 
40 113.6 3258.40 52.95 


表 3 描述 了 轨道 高 度 一 定 ,轨道 倾角 变化 对 应 
的 少 约束 转移 轨道 的 飞行 时 间 与 各 点 速度 增 量 的 变 
化 。 表 4 则 在 表 3 中 轨道 高 度 约 束 的 基础 上 增加 升 
交点 赤 经 约束 ,利用 分 层 微 分 修正 算法 ,分 析 了 轨道 
倾角 变化 对 应 的 多 约束 转移 轨道 的 各 项 参数 。 表 5 
给 出 了 固定 轨道 倾角 约束 , 升 交点 赤 经 变化 的 多 约 
束 转移 轨道 各 点 数值 。 由 表 3 与 表 4 可 知 , 对 于 固 
定 法 向 幅 值 Halo 轨道 的 固定 入 轨 点 , 随 着 轨道 倾角 
不 断 增 大 ,TTI 点 机 动 速度 增 量 变化 较 小 ,而 HOI 点 
施加 的 机 动 速度 增 量 变化 明显 。 表 4 与 表 5 表明 ， 
分 层 微 分 修正 算法 能 够 处 理 多 约束 问题 且 初 值 多 项 
式 能 够 较 好 地 保证 算法 收敛 。 
3.3” 初 值 多 项 式 对 不 同 幅 值 Halo 轨道 与 HOI 点 的 
适用 性 

在 目标 Halo 轨道 上 选取 不 同 HOI 点 对 飞行 时 
间 与 机 动 速 度 增 量 将 产生 影响 , 以 下 将 验证 本 文 所 
设计 的 初 值 多 项 式 对 于 不 同 HOI 点 在 转移 轨道 设 
计 中 的 有 效 性 。 

首先 将 Halo 轨道 按照 时 间 等 分 为 360 个 点 , 定 
义 Halo 轨道 与 会 合 系 xy 平面 相交 两 点 中 距离 地 球 


较 远 的 点 为 1 点 , 正 向 积分 一 周 后 对 Halo 轨道 上 点 


进行 编号 ,如 图 12 所 示 。 
Fery 90 


z/(10“AU 
y Ou 


-4 0.99 
0.989 x/((AU) 


12 目标 Halo 轨道 编号 定义 
Fig. 12 The definition of Halo trajectory numbering 


选取 1 点 作为 HOI 点 ,利用 初 值 多 项 式 设 计 了 
快速 转移 轨道 与 长 时 间 转 移 轨 道 。 快 速 转移 轨道 如 
13 所 示 ,此 转移 轨道 可 实现 天 文 任务 中 短 时 间 到 
达 目 标 轨道 的 目的 。 其 飞行 时 间 为 39.8 天 ,地 球 停 
泊 轨 道 逃 逸 速度 增 量 大 小 为 3262. 68 m/s, HOI 点 
机 动 增 量 大 小 为 224. 47 m/s。 图 14 则 描述 了 长 时 
间 条 件 下 的 平 动 点 转移 轨道 ,探测 器 施加 
3259. 03 my/s 从 地 球 逃 逸 后 , 飞行 185. 1 天 后 ,以 大 
小 为 58.78 m/s 的 速度 增 量 进入 目标 轨道 。 比 较 可 
知 , 对 于 1 点 入 轨 , 若 选择 短 时 间 飞 行 策略 ,需要 消 
耗 较 大 的 入 轨 速 度 增 量 , 而 长 时 间 飞 行 策略 则 将 增 
大 探测 器 飞行 时 间 。 


iD _ 上 上 An 


0 


JIO3AU) 


:9 
下 


-5 1 1 1 1 1 1 1 1 
0.988 0.99 0.992 0.994 0.996 0.998 1 1.002 1.004 
XAU 


13 ”快速 转移 轨道 
Fig. 13 The rapid transfer trajectory 


针对 不 同 HOI 点 ,当选 取 Halo 轨道 上 靠近 地 球 
一 侧 的 点 作为 入 轨 点 时 ,与 背离 地 球 一 侧 点 比较 , 虽 
然 需要 消耗 较 大 的 速度 增 量 ,但 能 够 缩短 探测 器 的 
飞行 时 间 。 因 此 ,本 文选 择 同 一 法 向 幅 值 的 Halo 轨 
道上 近 地 球 一 侧 点 作为 入 轨 点 ,所 设计 的 转移 轨道 
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AN 
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-0 
0.985 0.99 0.995 1 1.005 1.01 1.015 
vAU 


14 ”长 时 间 转 移 轨道 


Fig. 14 The transfer trajectory for long time 


结果 如 表 6 所 示 。 其 中 任务 约束 满足 : 停泊 轨道 高 
度 h = 200 km ,轨道 倾角 ; = 54° , 升 交 点 赤 经 
人 2=-15°*, 目 标 Halo 轨道 法 向 幅 值 4 = 
120000 km 。 从 表 6 可 见 ,针对 同一 幅 值 Halo 轨道 
不 同 HOI 点 ,满足 多 约束 条 件 的 转移 轨道 的 差异 主 
要 存在 于 HOI 点 速度 增 量 大 小 与 飞行 时 间 。 

表 6 不 同 HOI 点 转移 轨道 参数 


Table 6 Parameters of transfer trajectory for 


different HOI points 


HOI 点 “中 点 机 动 速度 。 HOI 点 机 动 速度 飞行 时 
” ， 坊 量 /(m *"s  ) ” 增 量 /(m*s7) 间 / 天 
140 3258.34 98.22 101.2 
160 3258.40 73.81 101.7 
180 3258.52 58.05 108.8 
200 3258. 63 54. 86 113.7 
220 3258.78 56.78 120.4 
240 3258. 89 67.24 126.9 


表 7 不 同 法 向 幅 值 Halo 轨道 设计 参数 
Table 7 Parameters of Halo trajectory design for different 


normal amplitudes 


法 向 幅 。 TTI 点 机 动 速度 ”HO0I 点 机 动 速度 飞行 时 
值 /km  ” 垃 量 /(m，s-!) 垃 量 /(m*s-!) 间 / 天 
50000 3216.61 66.86 109.8 
100000 3216.70 66.34 105.7 
120000 3216.77 65.83 104.6 
150000 3216.87 63.79 103.8 
200000 3216.98 58.11 102.9 
300000 3217.36 40.63 101.8 


表 7 选取 与 表 6 相同 的 地 球 停泊 轨道 约束 , 针 
对 不 同 法 向 幅 值 Halo 轨道 的 固定 入 轨 点 Node = 


180, 基 于 分 层 微 分 修正 算法 与 初 值 多 项 式 设计 了 满 
足 多 约束 条 件 的 转移 轨道 ,并 给 出 了 各 项 参数 的 数 
值 。 由 表 7 可 知 , 随 着 法 向 幅 值 不 断 增 大 ,转移 轨道 
飞行 时 间 与 HOI 点 机 动 增 量 呈现 减 小 趋势 。 

表 6 与 表 7 进一步 验证 了 本 文 设 计 的 初 值 多 项 
式 适 用 于 设计 不 同 幅 值 Halo 轨道 的 不 同 入 轨 点 的 
转移 轨道 。 


4 结 论 


本 文 针对 多 约束 条 件 下 Halo 轨道 的 转移 轨道 
进行 了 研究 ,利用 分 层 微 分 修正 算法 修正 转移 轨道 
末端 以 同时 满足 地 球 停泊 轨道 高 度 、 轨 道 倾角 、 升 交 
点 赤 经 及 航 迹 角 等 多 约束 条 件 。 并 且 , 对 迭代 改正 
量 与 状态 转移 矩阵 进行 了 变换 , 从 而 保证 与 加 快 了 
算法 收敛。 同时 ,针对 微分 修正 初 值 选取 与 敏感 性 
问题 ,通过 数据 拟 合 设计 了 能 够 为 控制 变量 提供 良 
好 初 值 的 多 项 式 。 初 值 多 项 式 能 够 满足 不 同 幅 值 
Halo 的 不 同 HOI 点 的 转移 轨道 设计 方案 。 通 过 具 
体 算 例 进一步 校 验 了 分 层 微分 修正 算法 在 处 理 多 约 
束 问 题 的 有 效 性 ,以 及 设计 的 初 值 多 项 式 的 实用 性 。 

本 文 在 仿真 过 程 中 虽然 取 定 了 特定 时 刻 的 地 球 
轨道 参数 ,但 对 于 其 他 时 刻 , 提 出 的 修正 算法 同样 适 
用 。 其 次 ,在 设计 转移 轨道 过 程 中 并 未 对 结果 进行 
优化 。 因 此 ,本 文 基于 分 层 微分 修正 算法 与 初 值 多 
项 式 设计 的 转移 轨道 是 满足 多 约束 条 件 、 转 移 时 间 
在 110 天 左右 且 相 比 于 短 时 间 转 移 , 能 量 比较 节省 
的 SEL1 点 转移 轨道 ,对 以 SELI 点 附近 Halo 轨道 为 
目标 轨道 进行 多 约束 转移 轨道 设计 时 具有 参考 
价值 。 
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